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Abstract 

To enhance the reflectivity of X-ray mirrors beyond the critical angle, multi- 
layer coatings are required. Interface imperfections in the multilayer growth 
process are known to cause non-specular scattering and degrade the mir- 
ror optical performance; therefore, it is important to predict the amount of 
X-ray scattering from the rough topography of the outer surface of the coat- 
ing, which can be directly measured, e.g., with an Atomic Force Microscope 
(AFM). This kind of characterization, combined with X-ray reflectivity mea- 
surements to assess the deep multilayer stack structure, can be used to model 
the layer roughening during the growth process via a well-known roughness 
evolution model. In this work. X-ray scattering measurements are performed 
and compared with simulations obtained from the modeled interfacial Power 
Spectral Densities (PSDs) and the modeled Crossed Spectral Densities for all 
the couples of interfaces. We already used this approach in a previous work 
for periodic multilayers; we now show how this method can be extended 
to graded multilayers. The upgraded code is validated for both periodic 
and graded multilayers, with a good accord between experimental data and 
model findings. Doing this, different kind of defects observed in AFM scans 
are included in the PSD analysis. The subsequent data-model comparison 
enables us to recognize them as surface contamination or interfacial defects 
that contribute to the X-ray scattering of the multilayer. 
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1. Introduction 

Multilayer coatings are known to enhance the reflectivity of extreme ul- 
traviolet (EUV), neutron, and X-ray mirrors at incidence angles larger than 
the critical one. As the energy of the incident beam increases, the smooth- 
ness of the surface becomes more and more important, because the roughness 
reduces the specular reflectivity and increases the X-ray scattering (XRS) in 
non-specular directions, leading to a degradation of the angular resolution. 
Depending on the speciflc application, periodic or graded multilayers are 
deposited using different techniques: anyway, to a variable extent, the depo- 
sition process triggers an evolution of the roughness from the substrate to 
the outermost layer. The interference of waves scattered at layer interfaces 
[l[ result in the final XRS pattern (Fig. [1]); therefore, in order to estimate 
the roughness impact on the Point Spread Function (PSF), a roughness mea- 
surement of all the multilayer interfaces would be needed. 

However, only the outer surface of the multilayer is accessible to direct 
topography measurements using, e.g., an Atomic Force Microscope (AFM). 
X-ray reflectivity (XRR) measurements as a function of the incidence angle, 
combined with a detailed flt routine to interpret the reflectivity scans (2!, Isl, 
allows a non-destructive, in-depth analysis of the multilayer stack structure 
(layer thickness in the stack, uniformity, smoothness), but does not enable 
the reconstruction of the PSD {Power Spectral Density) evolving throughout 
the stack. Nevertheless, an XRS computation based upon the sole thickness 
description and the outer surface PSD, assuming the rough topography to 
be exactly replicated in the stack, would in general return a diagram mis- 
matching experimental data (Fig. Hal) . 

In this paper, a modeling of the PSD evolution in the stack is used to 
compute the XRS diagram. A known multilayer growth model [3| provides 
the PSD growth across the stack, modeled from the measured PSD of the 
substrate and of the multilayer surface. The model physically describes the 
roughness of each interface as stemming from two effects in mutual compe- 
tition: the replication of the roughness of the underlying interface and the 
roughness introduced by the growth of the layer itself (Eq. [1]). As a result, 
the PSD increase from the substrate to the outer surface can be modeled by 
tuning the values of a few growth parameters ^ that can be tuned to flt 
the measured external PSD. Once the best-flt parameter values are set, the 
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"scattered": interfering X-rays at 9^ tB. 




i.j = 0, 1, ... 2N 



Figure 1: Scheme of X-ray scattering in a multilayer stack. Both "reflected" 
(i.e., in the direction specular to that of incidence) angle and "scattered" (i.e. 
in non-specular directions) rays result from the interference of elementary 
waves scattered at each boundary in the multilayer. 



internal PSDs for all the interfaces and the Crossed Spectral Densities (CSD) 
for all the couples of interfaces can be reconstructed. These quantities, in 
turn, determine the XRS expected from the multilayer. In particular, it is the 
cross- correlation between nearby or distant interfaces to affect the coherence 
between scattered waves, and, consequently, the amplitude of interferential 
features in the XRS diagram. The first-order perturbation theory (see (g], 0]) 
is used to compute the XRS diagram from the roughness PSD evolution in 
the stack gQ. 

Previous works js], have already implemented this model for periodic 
multilayers. However, broadband multilayers, like the ones in use in X-ray 
telescopes, have a graded structure ^] . In this work, we extend the formal- 
ism to graded multilayers and we apply it to two multilayer samples, a W/Si 
periodic sample and a Pt/C graded sample. Roughness analysis of the sub- 
strate and the outermost surface of the samples is performed with the Atomic 
Force Microscope (AFM) operated at INAF/OAB. The layer thickness mea- 
surement is obtained from the accurate fit jij of the XRR measurements 
performed with a BEDE-Dl diffractometer, also operated at INAF/OAB, at 
the X-ray energy of 8.045 keV (the Cu Kai fluorescence line). Eventually, 
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aiming at checking the correctness of the growth description, we have com- 
pared the expected XRS diagram to the one measured at with the BEDE-Dl 
at selected incidence angles, finding a very good agreement between the mod- 
eling and the experiment. Some preliminary results were already exposed in 
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a previous paper ^ ^ 

In Sect. [2] we retrieve the adopted growth model [3] and the XRS formal- 
ism applied to multilayers ^j. In Sect. Elwe describe the samples and the 
experimental setup. In Sect. IHwe show the modeling of the PSDs growth, 
as well as the predicted XRS vs. the experiments. For the periodic case, 
we have reanalyzed the data already treated [sj, showing that including the 
modeled PSD growth and adopting a more general electric field modeling 
leads to the best data/model matching. For the graded case, we show that 
the thickness trend that describes the stack and fits the XRS peak positions 
is univocally determined, and also matches the XRR measurement. Finally, 
we show how XRS can be a powerful tool to discriminate between surface and 
embedded defects, by including them in the PSD and checking if a proper 
XRS fit is obtained. The results are briefiy summarized in Sect. |3 A possible 
derivation of the formula used to modeljthe scattering diagram is sketched 
in 



Appendix A, or with more details in [9[. 



2. Modeling microroughness growth and X-ray scattering in mul- 
tilayers 

2.1. Microroughness growth model 

The roughness growth model solves a kinetic equation to describe the 
evolution of the rough profile z{x) with the thickness r of the film. For a 
single layer deposited onto a substrate, this equation reads 

^^-Hv»3Wk|1. (1) 

The model describes the roughening of the surface as a competition between 
a surface relaxation process and the increase in roughness due to the random 
nature of the deposition process. The relaxation process is parametrized 
with u and the positive integer n that varies with the kinetic mechanism 
that dominates the smoothing process 0]. The increase in roughness results 
from the deposition process and is described by a random shot noise term 77. 
The solution of Eq. [T]in terms of surface PSD [3] is 

p'-'V) = " '"ilr't:^'""' . (2) 

zz/|Z7rj I 
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where P™*(/) is the intrinsic bi- dimensional PSD of the layer surface, i.e., the 
PSD that the surface layer would have if the substrate were ideally smooth. 
This PSD is characterized by a plateau up to the maximum frequency corre- 
sponding to the cutoff wavelength /* = (z/r)", then decreases as a power-law 
of spectral index n. Vt represents the volume of the deposited atom, molecule, 
or nano crystal. 

When a stack of alternated layers is considered, the situation is compli- 
cated by the presence of two elements with different properties, i.e., different 
values of the parameters VL, u, and n. However, the formalism can be ex- 
tended by considering each single layer (whose upper surface is labelled with 
j = 0, 1,. . . moving from the substrate towards the surface) as growing 
upon its underlying layer, which acts as its substrate. In this way, one can 
write [3] -Po = Psubs and the PSD of the j^^ interface as a sum of the intrinsic 
contribution of the layer itself and of a term representing the rough profile 
partially inherited from the previous layer: 

PAf) = Pf^f) + P^^f) = P'T^f) + (3) 

In the second term of Eq. [3l 

a,(/) = exp(-H27r/rr,) (4) 

is a replication factor that describes to which extent the microrelief com- 
ponents of the (j — 1)*^ interface are replicated in the j*^ layer. Like the 
pint ^gj^j^i, the replication factor exhibits a cutoff at the spatial wavelength 
/* = (z/r)", that represents the transition above which surface features are 
damped out. At wavelengths larger than /*, the surface topography is al- 
most entirely replicated: the superimposition of the intrinsic term thereby 
triggers a progressive roughening of the multilayer interface, as a function of 
/ (Fig. |2]), to an extent depending on the values of the growth parameters. 
By fitting these values to the measured PSD growth from the substrate to 
the outer surface, the parameter values, and consequently the internal PSD 
evolution, can be determined. 

We have translated this formalism into an IDL-based program, firstly 
aiming at determining the fih, t'h, "z^-h, ^i, and n\ for the high- and low-Z 
density element, respectively. The values that best fit the growth from the 
substrate to the outer PSD, computed by recursive application of Eq. [3l are 
selected. Once the fitting parameters are found, we compute each Pj by 
iterating Eq. 12] j times from the substrate. The CSD between the generic 
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Figure 2: Transmission Electron Microscope section of a W/Si periodic mul- 
tilayer. The roughness profile is amplified throughout the stack ( Courtesy of 
C. Ferrari and L. Lazzarini, IMEM-CNR). 



j and the m interface, with j < m can be simply computed from the 
rephcation factors (Eq. H]): 

C,m{f) = aUf) ■ cim-i{f) ■ . . . ■ a,+i(/)P,(/) (5) 

The Pj's and the C^m's are the physical quantities that affect the intensity 
and the mutual coherence of the scattered waves, which in turn interfere to 
build up the XRS diagram. 

2.2. X-ray scattering 

To compute the XRS diagram at the X-ray wavelength A we apply the 
first order perturbation theory in grazing incidence, assuming the surface to 
be smooth and isotropic. For a single boundary characterized by a single 
PSD, the relation is a simple proportionality between the scattered intensity 
per angle unit and the mono dimensional PSD P{f) @, 0|, 

Tod/ " sm2^,sm^iP(/), (6) 

where the scattering angle 6s and the incidence angle d[, measured from the 
surface (Fig. [T]), are related to the X-ray wavelength A and to the spatial 
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frequency / as follows: 

i = , .. (7) 

/ I COS 6i — cos I 

The Qis polarization factor can be related to the reflectivity of the boundary 
i?F, computed via the Fresnel equations, 

Qis = y^MejMOi). (8) 

In the case of a stack of alternated layers, this relation is compli- 
cated by the interference effects among interfaces [H, @] . If the roughness is 
isotropic and the incidence angle is beyond the critical angles for total reflec- 
tion of both materials, the scattering diagram can be expressed in terms of 
the mono dimensional Pj's and C^m's using the following equation |9.], which 



represents a generalization of Eq. [6] (see Appendix A ): 

N 



_j=0 j<m 

(9) 

In Eq. [9] the layers are numbered from the substrate toward the surface. The 
spatial frequency / is still related to the scattering and the incidence angles 
via Eq. [71 a = 27r(sin^s -|- sin6'i)/A, Tj is the field amplitude transmittance 
in the j^^ layer, Ajm = (zj) — {zm) is the average distance between the j^^ 
and the m**^ interface, and 

XQ77-2 

^(A, e,, eo = -^Qis sin^ sin e, (lo) 

is the same proportionality factor appearing in Eq. |6l The polarization fac- 
tor is still provided by Eq. [HI being the intensity reflectivity of the single 
boundary between layers. Beyond the critical angles the values of Rp de- 
pend only a little on whether the reflection occur at the high-to-low density 
transition or the low-to-high one. 

It should be noticed that in Eqs. [9] and \1U\ both incidence and scattering 
angles are corrected for refraction in the stack. This is possible only if the 
incidence angle is beyond the critical one of both materials, otherwise the 
refraction angle will no longer be meaningful. Therefore, to maximize the 
intensity of the scattered beam, it is convenient to set the incidence angle at 
one of the first peaks after the critical angle of the multilayer. Secondly, it 
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should be kept in mind that Eq. M becomes inaccurate at scattering angles 
smaller then the critical one. 

An important point is the computation of the Tj coefficients. In a previous 
paper |5] we assumed an exponential decrease of the intensity throughout the 
stack [9l] when rays impinge at a Bragg angle. The method hereby adopted 
is applicable to both periodic and graded multilayers and makes use of the 
recursive theory of multilayer reflectivity to derive the transmittance of the 
partial stack of the outermost N — j layers. We have validated the recursive 



method by comparing its results with the findings of the IMD program [13 
finding a very good agreement. 

The result of Eq. [Hlwas added to a gaussian function at = 6i to simulate 
the specular reflected beam. Finally, the simulated XRS curve has been 
normalized to the flux impinging the sample, Jq, and smoothed to account 
for the finite detector acceptance angle dOg. 

The simulated scans are compared to the experimental XRS curves, ob- 
tained with the apparatus described in the next section. 



3. Experimental 

3.1. Samples 

In this work we consider two samples as test cases: 

1. a periodic multilayer with 40 bilayers of W/Si with dys/ = 19.3 A and 
dsi = 26.7 A deposited by e-beam evaporation onto a Silicon wafer 
(already characterized and analyzed 0]) and 

2. a graded multilayer with 100 bilayers of Pt/C deposited by magnetron 
sputtering onto a Silicon wafer, with layer thickness values as per the 
widespread supermirror design [ij], given by the power law d{j) = 
a{b + where d is the d-spacing and j is the bilayer index from the 
top of the stack. The nominal power-law parameter values are 

Pt : a = 31.0 A,b = -0.94, c = 0.23, 
C: a = 53.0 A, b = -0.88, c = 0.21 



being Platinum the ffist layer deposited onto the substrate. The out- 
ermost layer is Carbon. 
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3.2. Atomic Force Microscope 

The surface topography of the samples before and after coating was mea- 
sured with the stand-alone Atomic Force Microscope (AFM) operated at 
INAF/OAB, a Veeco Instruments, mod. Explorer®. The scanner can image 
surface areas of 1/xm to 100 /im side. The AFM was operated in tapping 
mode, using an Antimony- doped Silicon probe with a proper frequency in 
the 100-250 kHz range. The maximum lateral resolution is 4 nm and the 
vertical resolution is better then 1 A. 

Several scans were performed on each sample to check the uniformity of 
the deposition process and the possible presence of different kinds of defects. 
From 10 /im and 1 fim wide scans we computed the PSDs of the substrate 
and the outer surface of the multilayer in the spatial wavelength range 8 nm 
-10 /im. This spectral range is also typically involved in the roughening 
process [sl]: at larger wavelengths, the rough profiles are essentially replicated 
without relevant roughness increase, whilst high frequencies are damped out. 

3.3. X-ray Reflectivity 

The thickness values of the multilayers under test were investigated via 
X-ray Reflectivity (XRR) measurements with the BEDE-Dl diffractometer 
from Bede Scientific Instruments Ltd®. A Cu anode X-ray tube is used as 
source. A channel-cut Si crystal monochromator and a slit collimator are 
used to filter and shape the Cu Kal fluorescence line at 8.045 keV: the first 
slit (50 nm wide) removes the X-ray continuum and the Ka2 line diffracted 
by the monochromator, whereas the second slit (10 fim wide, very close to 
the sample) reduces the beam width for the sample to collect it entirely, 
also at very small incidence angles. The reflected beam is then collected 
through a 800 /xm wide slit by a scintillation detector with a pulse-height 
discriminator that limits the intrinsic background to < 0.2 counts/sec. The 
measurement is performed via a 6 — 26 scan with 20 arcsec steps, on the order 
of the achieved beam divergence (15 arcsec). After the scan, the reflectance 
is normalized to the incident flux. The detailed fit of the XRR curve, to 
determine the stack parameters (thickness, density, roughness), is performed 
using the PPM program 

3.4. X-ray Scattering 

The BEDE-Dl diffractometer was also used to measure the scattering of 
8.045 keV (1.541 A) X-rays at incidence angles corresponding to the first 
peaks after the critical angle. In this double channel-cut Si crystal 
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Scintillation 
detector 



Figure 3: Scheme of the XRS experimental equipment. 

monochromator was used in order to improve the beam coUimation (to within 
10 arcsec FWHM). In this configuration, a thin slit is not needed to isolate the 
Cu Kal line. On the other hand, the incoming fiux has to be intense in order 
to scan the scattered beam up to angles as large as 15000 arcsec. Moreover, 
the incidence angle will be near 1 deg from the surface; therefore, the beam 
does not need to be very thin to be entirely collected on the sample, so a 
400 /im-wide slit (wider than the one used for the XRR) was positioned after 
the monochromator. A 400 fim slit in front of the detector, at a 340 mm 
distance from the sample, provides a sufficiently small angular acceptance 
for our scopes {d6s ~ 320 arcsec FWHM). The measurement of the scattered 
intensity, dig (Eq. IH]), was performed by scanning the detector in the plane 
of incidence with 20 arcsec steps. The layout of the experimental apparatus 
is sketched in Fig. |3l 

Application of Eq. [7] shows that a detector scan with 6s G [~ 3000 -v- 
15000] arcsec, with 9i = 2200 deg, at A = 1.541 A, corresponds to spatial 
wavelengths between a few microns and fractions of micron, i.e., the ones in 
the sensitivity window of the AFM (Sect. 13.21) and mostly involved in the 
roughness growth process. 
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4. Results vs. modeling 

4.1- Sample A 




Scattering angle (arcsec) Scattering angle (arcsec) 



(a) (b) 

Figure 4: Experimental (red lines) X-ray scattering for the sample A vs. 
modeling (black lines), (a) Assuming a constant PSD in the stack returns 
a poor modeling, (b) Accounting for the roughness evolution through the 
stack reproduces more correctly the experimental XRS curve. 

The complete characterization of this sample was already reported in a 
previous work jsj. The roughness growth was also studied in detail therein. 
The XRS scan, taken at the 1*^* Bragg angle of incidence, exhibits an apparent 
peak at ~ 10000 arcsec, stemming from the constructive interference of 
scattered waves. In |5|], however, the computation was performed using an 
exponential trend for the Tj coefficients (Sect. I?]2|) . which can be only applied 
to periodic multilayers in Bragg incidence. We have reconsidered those data 
to validate the general method used heretofore to compute the Tj's. The 
experimental XRS scan vs. the result of the modeling out of Eq. IHlis displayed 
in Fig. m In order to show the impact of the PSD evolution in the XRS 
diagram, we have preliminarily assumed no evolution of the PSD throughout 
the stack with a complete correlation {Cjm = Pj = -Pouter) at all frequencies. 
The modeling clearly overestimates the measurement (Fig. Hall : as expected, 
this denotes a gradual evolution of the roughness from the substrate from 
the outer surface. Moreover, the hypothesis of a complete correlation at all 
frequencies {aj = 1 for all j, Eq. H]) endows the simulation with a second 
XRS peak near 17000 arcsec, which is not observed in the experimental XRS 
curve. In contrast, setting the Spectral Density functions resulting from the 
correct modeling of the roughness growth yields a modeling in much better 
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accord with the experimental scan (Fig. l4b|) . The model-data matching also 
proves the correct trend of the electric field coefficients throughout the stack, 
computed using the general method. 

4-2. Sample B 

For this graded multilayer sample we have directly measured with the 
AFM (Sect. 13. 2p the roughness of the outermost Carbon layer in the 10 /xm - 
5 nm spatial wavelength range, while the substrate roughness was supposed 
to be the same of a standard Silicon wafer [l^, e.g., 2 A rms in that spec- 
tral window. The multilayer surface, as measured with the AFM, exhibits 
crowded point-like defects in ejection (Fig. [Sa|) that increase the rms of the 
external surface to 3.7 A in the same spectral range. Besides, other defects of 
bigger size increase the rms from 3.7 A to 7.3 A in the 10 /im scan (Fig. [5b|) . 
These bigger defects are homogeneously distributed over the sample, but it 
is difficult to ascertain if they are surface contaminations or they stem from 
the roughness growth itself. In the latter case, such defects would be effec- 
tive for X-ray scattering, whilst in the former one, 8.045 keV X-rays would 
be almost unaffected by their presence. For this reason, we have performed 
two modelings of the roughness growth, including or excluding the biggest 
defects. The two PSD evolutions are shown in Fig. [5c] and [5dl Both models 
were used to simulate the expected XRS diagram: the comparison to the 
experimental curve will show that the modeling of Fig. |5d]is correct. 

Prior to the scattering measurement, the XRR of the sample at 8.045 keV 
has been measured (Sect. 13.31) . The PPM program was subsequently apphed 
to perform a detailed fit of the refiectivity 01 and return the best-fit power- 
law parameters for the stack (Tab. [H fit 'A'). The layer thickness values 
were used to firstly model the roughness growth (Eq. [3]) and to subsequently 
simulate the XRS diagram (Eq. ^ of X-rays impinging at the second peak 
after the critical angle (^i ~ 3000 arcsec). The electric fields were computed 
using the recursive method, which is the only one applicable for a graded 
multilayer and already tested with the Sample A (Sect. 14. ip . We already 
note that adopting the PSD evolution modeling that includes all the sur- 
face defects (Fig. [5dl) yields the correct scattering level. Nevertheless, the 
XRS diagram simulated in this way did not perfectly fit the measured XRS 
(Fig. [Taj); for instance, the peak positions did not match. This denotes some 
departure of the actual thickness trend from the modeled one. 

We have so found that a slightly different power law (Tab. [H fitB) matches 
much better the scattering peak positions (Fig. [7b|) . This new power- law 
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(c) (d) 

Figure 5: Roughness analysis from AFM scans of the outer surface of the 
sample B. The AFM image (a) shows that the surface is crowded with point- 
like defects in ejection, and also exhibits some defects of bigger size, (b) 
The bigger defects can be included or excluded from the PSD computation. 
The PSD evolution in the stack was consequently computed (c) with and (d) 
without the major defects. The XRS measurement shows (Fig. [8b|) that the 
correct modeling is (d). 
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Figure 6: XRR angular scan at 8.045 keV of sample B, compared to the 
reflectivity computed from two slightly different stack models. They both 
match the XRR measurement, but only the fit 'B' is also in agreement with 
the peak positions seen in the XRS scans (Fig. [Tbll . 




Figure 7: Experimental XRS of the sample B (red lines) at 9\ = 3000 arcsec, 
compared with the simulated ones from the PSD evolution, including the 
major defects (Fig. [5dl black lines), (a) XRR fit 'A', (b) XRR fit 'B'. The 
layer thickness values provided by the fit 'B' of Tab. [Hreproduces much better 
the observed peak positions. 
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actually returns layer thickness values that do not differ by more than 3 
A from the first one, but in the outermost C layer, which turns out to be 
15 A thicker. Quite surprisingly, the XRR measurements actually matches 
both A and B power laws (Fig. |6]). A possible interpretation could be that 
the specular reflectivity at 8.045 keV is a little sensitive to the thickness of 
the outer C layer, because most of the XRR curve features that drive the 
fit program are determined more by the multilayer internal structure, rather 
than by the thickness of the capping layer. In contrast, an XRS measurement 
can be more sensitive to the increased roughness generated by a thicker C 
layer. The experimental scattering diagram is now well reproduced by the 
modeling {xl ~ 0.15). We therefore assumed the power law 'B' of Tab. [T]to 
be the correct one. 



Table 1: Fitting power law parameters 



Case 


Pt: a 


[A], b, c 


C: a 


[A], b, c 


Fit A 


31.08, - 


-0.90, 0.30 


53.16, 


-0.96, 0.19 


Fit B 


31.00, - 


-0.94, 0.23 


53.00, 


-0.88, 0.21 



Modeling the stack with the parameters given by the fit B, we hereafter 
consider the effects of including different defects in the PSD, and therefore 
in the XRS computation. If the sole point-like defects are included (Fig. I5a|) 
and the growth parameters are tuned accordingly, the peak heights are un- 
derestimated in the simulated XRS diagram (Fig. I8a|) . Conversely, if all the 
visual defects in the AFM maps are used to compute the final PSD of the 
multilayer, the scattering diagram matches much better the experimental 
XRS curve (Fig. [8b|) . We therefore conclude that the visual defects other 
than the point-like defects observed on the multilayer surface (Fig. [Sal) are 
related to the deposition process, and not just dust contaminations. 

A further validation has been obtained from an XRS measurement at a 
different incidence angle, i.e., the first XRR peak after the critical angle (^i = 
2200 arcsec). The comparison between theory and experiment (Fig. [9]) shows 
a good accord also in this case {xf, ~ 0.29), even if the peak positions are not 
exactly reproduced: the stack structure might be slightly inhomogeneous 
and the two XRS measurements, performed in two different runs, might 
have been performed at different locations. Consequently, also the power 
law parameters might require some adjustment. 
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Figure 8: Experimental XRS for sample B (red lines) at Q\ = 3000 arcsec, 
assuming the fit 'B' to model the stack structure, compared with the sim- 
ulated ones from the PSD evolution (black lines): (a) considering only the 
point-like defects of Fig. |5al (b) considering all the visual defects in the AFM 
maps. 
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Figure 9: Experimental XRS for sample B (red lines) at $[ = 2200 
suming the fit 'B' to model the stack structure, compared with the simulated 
ones from the PSD evolution (black lines) including all the visual defects in 
the AFM maps. 
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5. Conclusions 

The XRS sensitivity of a multilayer to the roughness evolution and to the 
actual thickness trend in the stack makes it a powerful diagnostic instrument, 
provided that it is coupled to an opportune modeling of the PSD evolution 
and that the XRS pattern can be reliably modeled (Eq. 4). So far js], [l^ 
the modelization of XRS for multilayers was limited to the periodic case. 
In this paper we have shown that the formalism can be extended to graded 
multilayers. In the tested cases we have reached a good experiment-modeling 
simultaneous matching for both PSD growth and X-ray scattering. The com- 
parison between XRR and XRS fitting shows that for graded multilayers the 
XRS is more sensitive than XRR to the actual thickness trend in the stack, 
since even a few angstroms variation significantly changes the XRS peak po- 
sitions. Moreover, by fitting the experimental XRS, it becomes possible to 
discriminate defects that produce scattering from surface contaminations, an 
analysis impossible to obtain from the sole AFM scans. 
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Appendix A. X-ray scattering in a multilayer coating 

In this appendix we briefly derive Eq. Oused in Sect. l2.2l to describe the 
X-ray scattering from a multilayer. The complete computation can be found 
in [9j. Let the multilayer substrate surface to be a rectangle of sides Li and 
L2 in the xy plane (Fig. lA.lOp . with the coating that covers a surface S at 



z = and occupies a volume V at z > 0. The electric field of amplitude Eq 
impinges the sample in the xz plane at a shallow incidence angle 6i and partly 
travels across the coating. The electric field E{z) amplitude varies with the 
depth in the stack owing to the progressive reflection and absorption in the 
multilayer. The scattered intensity results from the interference of elementary 
waves scattered by electrons in the entire stack. 

We denote with N'lz) the number of free electrons per volume unit, as a 
function of the depth in the coating. This number is related to the material 
density p, the atomic weight A, the scattering coefficient /i, as J\f{z) = 
pN^fi/A, where N\ is the Avogadro number. The scattered power dig into 
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Figure A. 10: Reference frame used in the computation of the XRS diagram 
of a multilayer coating. 



the solid angle dfls 
expressed as: 

dL 



cos^sd^sd0s at the polar coordinates 



c da 



can be 



(A.l) 



where c is the speed of light, is wave vector of the incident wave, k is 
wave vector after the scattering, r = (x, y, z) is the location of the scattering 
electron, and da /dVts = ^r'^{l + cos^(A6')) is the Thomson cross section for 
unpolarized radiation, with A6 the angular deviation and = 2.8 x 10~^^ cm, 
the classical electron radius. Substituting the expressions of the wave vectors, 
we define 



fx 

fy 
a 



cos 9s cos 0s — cos 9i 

A ' 

cos sin 0s 



27r 



A 

sin ^s + sin 6i 
A ' 



(A.2) 
(A.3) 
(A.4) 
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where the angles are evaluated at a proper average value of the refraction 
angles inside the stack. Using Eqs. IA.2I through IA.4t we can rewrite Eq. lA.ll 
as follows: 



d/s 



c da 



^2^g-2«/.x / E{z) M{z) e 



(A.5) 



where x = (x, y) and / = {fx, fy). We now define zq, zi, Z2, ■ ■ ■ ,zn to be the 
boundary profiles between the N layers from the substrate {zq = 0) to the 
outer surface. We also denote with (2:0), {zi), {Z2), ■ ■ ■ , {zjy) the average profile 
heights. Since the multilayer materials are supposed to be homogeneous, the 
electron density is constant in each layer and can take on the values A/h 
or A/i, for the high-density or the low-density material, respectively. If the 
absorption is negligible, also the electric field amplitude is nearly constant in 
each layer; therefore the integral in z can be solved. 



/2]V 1 ^ 

dzE{z)^{z) e"'"^ = —J2 {^j+iEj+ie 
j=0 



- MjEjC-''"'^) . (A.6) 



Substituting this result into Eq. IA.5t and re-arranging the terms we obtain 

2 



dJs 



c dcr 



N 
j=0 



la 



d^x e"^'^*^- e~*"^^' 



s 



(A.7) 



where Mn+i ~ 0, i.e., the electron density of the ambient. In most cases 
the contribution of the uppermost interface is negligible, so the difference of 
electron density at the alternated interfaces can be written as (— l)-'(A/'h— A/i). 

Owing to the smooth surface approximation, the exponentials can be 
approximated as e~*"^^ ^ e~*"^^^^(l — ia/S.Zj\ where ^Zj = Zj — (zj). Sub- 
stituting and developing the computation j9[, the O*'^ order term yields a 
delta-like term that represents the ray reflected at / = 0, i.e., in the specular 
direction (0s = 0, 9s = Oi). The subsequent non-zero term describes the 
scattering at the first order: 



1 dJs _ (ATh - MY da 
lodQs L1L2 sin6'i dfis 



N 
j=0 



-ia(zj) 



d^x Azj e 



(A.8) 



where Tj = Ej/Eq and the diagram was normalized to the incident intensity, 
Jo = {c/8tt)LiL2Eq sinOi. The remaining surface integrals are the Fourier 
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transform of the rough profiles of the layer boundaries, then executing the 
products the expression becomes 

1 dL (AAfV da 



Iq dflg sin 9i dflg 



' N 

_j=0 j^m 

(A.9) 

In Eq. IA.9I we have defined Aj^ = \ (zj) — (zm) \ , A.J\f = Ah — A/i- Pj and Cjm 
are the 2D Power Spectral Densities and the Crossed Spectral Densities in 
the stack, function of and fy. 

If the roughness is isotropic, in grazing incidence the scattering diagram 
essentially lies in the incidence plane 0]: so we integrate Eq. I A. 91 over (ps- 
Since the scattering power is essentially concentrated near 0s ^ 0, the vari- 
ation of the cross-section with 0s is negligible. From Eqs. IA.2I and IA.3I we 
have cos 9s d(/)s ~ A dfy if (ps ~ 0. So the integration operates on the sole 
Power Density functions, and we have 

+ (A,10) 

Jo dt^s sm Ui di Is 
where the cross section is evaluated at (^^s + ^i) and we have set 

N 

W^) = H^'PsU.) (A.ll) 

j=0 

P2{U) = 2 5^(-l)^+'"T,T^cos(«A,„,)C,„(/,). (A.12) 

Using now the expression of the optical constant 5 ^], we can express 
AM in terms of the change of 5 at the interfaces, A5. Accounting for the 
angular dependence of the electronic cross section, the proportionality factor 
in Eq. lA.lOl becomes 

It is easily checked, e.g. via a series development, that 

1 + cos2(^i + Os) '^S'' (1 + cos^ 2e,fl\l + cos^ 2esf'\ (A.14) 

Finally, using Eq. IA.14[ and recalling the approximated form of the Fresnel 
equation ^ for the grazing incidence reflectivity E^'^iO) = A5/(2sin^6') at 
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the generic angle 6 beyond the critical one, the proportionality factor turns 
into 

1 r 2 

K{\,9,,9,) ^ ^sine.sm' e,[RF{9,)RFi9s)]'^' . (A.15) 

By substituting the expression of Eq. IA.15I into Eq. lA.lOl and dropping the 
subscript 'x' in the spatial frequency because of the sample isotropy we even- 
tually obtain Eq. |9l 
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